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Q\ • Abstract 

1^ I We discuss a recent analysis of tt threshold effects and its implications for the 

I ' determination of electroweak parameters. We show that the new formulation, when 

Oh! 



applied to the p parameter for mj ~ 158 GeV, gives a result of similar magnitude to 
(~| , those previously obtained. In fact, it is quite close to our "resonance" calculation. 

We also present a simple estimate of the size of the threshold effects based on an 
elementary Bohr-atom model of ti resonances. 



1 Introduction 

The analysis of threshold effects involving heavy quarks and their contribution to the 
determination of electroweak parameters has been the subject of a number of studies in 
the past |1|J^,|,||,^. Very recently, F.J. Yndurain has discussed, in the framework of the 
Coulombic approximation, ti threshold effects on the renormalized vacuum-polarization 
function n(s) — 11(0) associated with conserved vector currents [^]. Here n(s) is the un- 
renormalized function defined according to n^j,(g) = {q"^ 9 fiu — (l^i(lu)^{(l'^) , with s = q^. For 
rrit = ^/SiTLz ~ 158 GeV and s = m|, the author of Ref. finds that the threshold effects 
are significantly smaller than the perturbative 0{as) calculation. From this observation 
he concludes that threshold effects are generally small for s <C 4mj and that the "large" 
results reported in Ref. |Q concerning their contribution to electroweak parameters are 
not supported by his "detailed, rigorous calculation." However, ti threshold effects influ- 
ence electroweak parameters chiefly through the p parameter. The function H(s) — H(0), 
although very important in its own right, has very little to do with the effects discussed 
in Ref. f^. In fact, it is well known that heavy particles of mass m'^ ^ s decouple in 
this amplitude. For this reason, it should be obvious that the leading effects discussed in 



Ref. 1^ do not arise from this amplitude. Thus, conclusions drawn on the work of Ref. |^ 
from the discussion of n(s) — n(0) are without foundation. In order to show in the sim- 
plest possible way what the correct conclusions ought to be, in Section 2 we apply the 
formulation of Ref. [^ to the study of leading threshold contributions to the p parameter 
for rrit ~ 158 GeV and other values of m^. In principle, this requires only minor modifica- 
tions of the relevant formulae of Ref. 0, mainly in the prefactors. However, for reasons 
explained in that Section, we find it necessary to re-evaluate the threshold corrections of 
that paper. We then find that, contrary to the conclusions of Ref. Q, this leads to results 
similar in magnitude to those reported in Ref. p. In fact, the approach of Ref. [§ gives 
threshold corrections to the p parameter quite close to our own "resonance" calculation. 
In order to make more transparent the size of these leading threshold corrections relative 
to the perturbative 0{as) calculations, in Section 3 we present a simple estimate based 
on an elementary Bohr-atom model of toponia. We also briefly comment on the shifts 
induced in electroweak parameters and the magnitude of 0{a1) corrections. 

2 Threshold corrections to the p parameter 

For large irtt, the dominant threshold effects discussed in Ref. can be related to cor- 
rections to the p parameter. As the current iptl^'^t is conserved and the p parameter is 
defined at s = 0, it is clear that tt threshold effects arise in this case from the contributions 
of the axial- vector current. As explained in Ref. |^, on account of the Ward identities 
such contributions involve '^m X^{s,mt,mt), where A"^ is a longitudinal part of the axial- 
vector polarization tensor 11^^. Furthermore, for non-relativistic, spin-independent QCD 
potentials, — Qm \^ can be identified to good approximation with Q^mll (in Ref. 0, 11 
is called Il^{s)/s). In particular, we note that both amplitudes receive contributions 
from nS states. One then finds that, in the formulation Ref. [^, the leading threshold 
correction to the p parameter is given by (cf. Eqs. (5.1, 5.2a) of Ref. P) 

SiAp)ti,r = --^ fds' Qmlithris'). (1) 

Here Q'mHt/ij.(s') denotes contributions from the threshold region not taken into account 
in the usual perturbative C(as) calculation. 

Our aim is to employ the analysis of Ref. to calculate Eq. (^ and then to compare 
the answer with our results. In the formulation of that paper, based on the Coulombic 
approximation, there are two contributions to ^mJ\thr[^s'\ One of them arises from 
the toponium resonances and is called 3?7iHpoie(s'). The other represents a summation 
of (a^/f )" in = 1,2,.. .) terms, integrated over a small range above threshold, after 
subtracting corresponding 0{as) contributions. The first one can be obtained from the 
expression 
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where N^ = 3, Cp = {N^ - l)/(2A/'c) = 4/3, n/ = 5, /?o = H - 2n//3 = 23/3, 
M„ is the mass of the nS toponium resonance in a Coulombic potential 
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Cpal{li)m^/{2n^) is the square of its radial wave function at the origin, and a;s(yu) = 
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(/i)(l + ba,ifi)/7i), with b = 7i?(llA^c - 2n/)/6 + (31Ar, 
1^. Inserting Eq. (^ into Eq. (|1]), we obtain 
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where Xt = {N ^G pml / Stx"^ \/2) is the Veltman correction to the p parameter [^. In the 
numerical evaluation of the threshold corrections to n(s) — 11(0) reported in Ref. P], 
/i is chosen to be m^, rrit = \pimz is assumed, and otsivnz) = 0.115 ± 0.01 is taken. 
Unless stated otherwise, we shall adopt these values in the following. Furthermore, 
{15/16)7i({3)Cpa^{fi) is found to equal 1.81 x 10~^, the summation J2'^=2 i^ neglected, and 
{1 + {3Poas / 2TT)[\iii{fi/ C Fdsfnt) — 1e]} is given as 1.08. It is apparent that the numerical 
value given in Ref. for the last factor is too low. The correct value is 1.315. Inclusion 
of the neglected sum raises the expression between square brackets in Eq. (^ to 1.437, 
instead of 1.08 0. This leads to 
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The contribution from the small range above threshold can be gleaned from Ref. ^. 
Using V = {l — Amf/sY^'^ as integration variable {v is the top-quark velocity in the center- 
of-mass frame) and approximating (1 — v"^)^"^ ^ 1 in the integrand, the contribution to 
Ap equals —{16/15)6thrXt, where 6thr is a quantity studied in Ref. 0. For s <^ 4mj, 6thr 
can be written as 
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where an overall factor (l—v'^/3) has been omitted under the integral. In Eq. (^), the term 
involving B{v) represents a summation of (vrCi^as/f )" contributions valid for large values 
of this parameter, while the last two correspond to the subtraction of the perturbative 
calculation up to 0{as) in the small-w limit. We have evaluated the latter at rrit, as this 
is demonstrably the proper scale to be employed in the perturbative calculation 0. We 
note that the integrand of Eq. (^ is renormalization-group invariant through Ola"^) |T0| . 
In Ref. 0, the value vq = 7rCFOis{m't)/V^ ~ 0.314 is chosen, the exponential and terms 
of higher order in as are neglected, and the answer given as 
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Although the analytic suinination of {ttCfo^s/v)"' terms is theoretically interesting, there 
are unfortunately a number of problems in the evaluation of 6thr carried out in Ref . P| : 



1. Evaluation of Eq. (0) as it stands gives a negative result, 6thr = —3.82 x 10 '^ jTl 



which obviously contradicts the well-known fact that multi-gluon exchanges lead to 
an enhancement of the ti excitation curve. 

2. Equation (|^ is not renormalization-group invariant, as the coefficient of ^^(/i) does 
not match correctly that of In/i. Subject to the approximations explained before, 
the correct, renormalization-group invariant expression is obtained by including an 
additional term [as(/i) — as{mt)]/2 within the parentheses of Eq. (|^. We note in 
passing that this additional term may be traced to a change of scale in the last term 
of Eq. d^), from the arbitrary value /x to the proper physical choice rrit. For fi = niz, 
the value of the corrected expression is 6thr = —3.83 x 10""^, i.e., essentially zero, 
and very different from the value reported in Ref. [^]. 

3. As we have a near cancellation of relatively large contributions and, moreover, the 
result should be positive, it is clear that, for the chosen value of vo, the neglect of 
the exponential in Eq. (||) is not justified. In fact, evaluating numerically Eq. (H) 
with the exponential included, we find 6thr = 1-55 x 10~^. This value is positive, as 
it should be, and much larger than the answer obtained without the exponential. 
Actually, it is 2.9 times larger than the result reported in Ref. 0. 

We also note that the integrand of Eq. (H) vanishes at 1.002 t>o, where Vq is defined above 
Eq. (0). Thus, Vq is a reasonable value to use as the upper limit of integration because at 
that point the resummed series and the perturbative 0{as) contributions nearly coincide; 
evaluation of the integrals up to the point where the integrand actually vanishes leads 
to negligible changes. On the other hand, a possible weakness of the method is that vq 
is rather large and for such values of v it is not clear that the resummed expression is 
valid. Nevertheless, for our present purpose, which is the evaluation of the contribution to 
Ap according to the prescriptions of Ref. 0, we use the above value of Vq. For the same 
reason, we use fi = rriz, although this is not a characteristic scale in connection with the p 
parameter; using the appropriate value, p = rrit, would make the radiative correction and 
the overall result slightly smaller. However, we include the effect of an additional overall 
factor (1 — f^/3)/(l — f ^)^, which should be appended to the integrand of Eq. (j^) in the 
case of Ap; it increases the result by about 4.3%. We then find that the contribution to Ap 
from the range < f < iiq above threshold is —(16/15)1.55 x 10"^ ■ 1.043 x^ = — 0.0172 x^. 
Combining this result with Eq. (P, we find that, after correcting the errors discussed 
above, the formulation of Ref. leads to S{Ap) = —0.0450 Xt- For comparative purposes, 
we rescale this result to the case as{mz) = 0.118, the value used in our calculations, and 
obtain, for rrit = VStuz, 
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5{Ap)tf,r = -0M86xt (Ref. [6]). (8) 

In our own work we have applied two different methods to evaluate the imaginary parts 
near threshold. The first one is the resonance approach of Ref. 0, which assumes the 
existence of narrow, discrete tt bound states characterized by Rn{0) and M„. Moreover, in 



Ref. [0 a specific interpolation procedure is developed to implement the matching of the 
higher resonances and the continuum evaluated perturbatively to 0{as)- The second one 
is the Green-function (G.F.) approach of Ref. 0, which takes into account the smearing of 
the resonances by the weak decay of its constituents and leads to a continuous excitation 
curve. Both approaches make use of realistic QCD potentials, the Richardson and the Igi- 
Ono potentials, which reproduce accurately charmonium and bottonium spectroscopy and 
are expected to describe well toponia, too. These QCD potentials contain a term linear in 
the inter-quark distance, r, to account for the confinement of color. Detailed studies reveal 
that the shape of the Green function is not very sensitive to the long-distance behaviour 
of the potential [|1^. This may be understood by observing that the top quarks decay 



before they are able to reach large distances. The rapid weak decay of the top quarks, 
which causes the screening of the long-distance effects, is properly taken into account 
in the Green- function approach of Ref. 0], while it is not implemented in the resonance 
approaches of Refs. |^ and [Q. In the latter case, it is clearly more consistent theoretically 
and more realistic phenomenologically to keep the linear term of the potential, as is done 
in Ref. ^ but not in Ref. |^. On the other hand, both resonance and Green-function 
approaches of Ref. effectively resum the contributions of soft multi-gluon exchanges in 
the ladder approximation |]3|JT3[|. This automatically includes the final-state interactions 
emphasized in Ref. 0. However, we stress that all these methods are based on a non- 
relativistic approximation. There are additional contributions due to the exchange of hard 
gluons, which give rise to sizeable reduction factors [l^, e.g., (1 — SCpas/ii) in the case 

ofApi- 

Because of the more complicated potentials used in our two approaches, we have to 
rely on numerical computations. For m^ = ^/Smz and as{mz) = 0.118, we find 

S{Ap)thr = -(0.034 ±0.010)xt (G.F.), (9) 

SiAp)thr = -(0.042 ±0.013)xt (res.), (10) 

where we have included the 30% error estimate given in Ref. [^. The corresponding 
perturbative 0{as) contribution [|l5| is, for rrit = y/STUz, 



5iAp).. = — ^^(y + ij-* 

= -0.0991 Xi. (11) 

Unlike 6{Ap)thr, 6{Ap)a^ obtains important contributions arising from the non-conserved 
vector and axial-vector currents associated with the PF-boson vacuum-polarization func- 
tion. 

It is apparent that the result for 5{Ap)thr obtained in the formulation of Ref. |^ 
(Eq. (H)) is of the same magnitude as our two evaluations (Eqs. (^, (|l3))- It amounts 
to 49% of the 5{Ap)oi^ correction, while our results of Eqs. (^, (jT^) correspond to 34% 
and 42%, respectively. Thus, it is somewhat larger than our resonance calculation and 
significantly larger than our G.F. result. Part of the difference is due to the fact that 
we have included a hard-gluon correction ||1^, (1 — 3Ci7'as(M„)/7r) (cf. Eqs. (4.2b, 4.3d) 



of Ref. p), which has not been incorporated into Eqs. (H), (§). Note that Eq. (^ must 
not be multiplied by this factor, since only non-relativistic terms are subtracted in the 



TABLE I. ti threshold effects on Ap relative to the Veltnian correction, —100 x 
6{Ap)thr/xt, as a function of m^. The calculation based on Ref. p (total [6]) is compared 
with our previous resonance (res. [5]) and G.F. (G.F. [5]) results. For completeness, the 
contributions from below (pole [6]) and above (thr. [6]) threshold are also displayed sep- 
arately in the first case. The hard-gluon correction is included in the three calculations 
and the input value as{mz) = 0.118 is used. 



mt [GeV] 


pole [6] 


thr. [6] 


total [6] 


res. [5] 


G.F. [5] 


120.0 


2.84 


1.99 


4.83 


4.81 


3.58 


140.0 


2.73 


1.92 


4.65 


4.43 


3.43 


157.9 


2.64 


1.86 


4.50 


4.17 


3.36 


180.0 


2.55 


1.79 


4.35 


3.91 


3.35 


200.0 


2.47 


1.75 


4.22 


3.71 


3.41 


220.0 


2.41 


1.70 


4.11 


3.53 


3.51 



TABLE II. Perturbative 0{as) and ti threshold |^ contributions to Ap relative to the 
Veltman correction, —100 x S{Ap)aJxt and —100 x S{Ap)thr/xt, for nit = y/Smz as a 
function of Ppert- The input value as(m^) = 0.118 is used. 



f^pert 


pert. 0{as) 


total [6] 


sum 


mt/2 
mt 
2mt 


10.96 
9.91 

8.87 


4.05 
4.50 
5.18 


15.00 

14.42 
14.05 



integrand of that equation. If this correction is applied to Eqs. (0), (|[), 6{Ap)poie becomes 
—0.0245 Xt and, rescaled to as{mz) = 0.118, the overall result in the formulation of Ref. 
is — 0.0450 a;t, instead of Eq. (||). This is quite close to our resonance calculation (Eq. |TD|). 
For the reasons explained in Ref. (see also the discussion), for values of rrit > 130 GeV 
we have expressed a preference for the G.F. approach. On the other hand, the three 
calculations amount to only 3.4%-4.5% of the 0{a) contribution, Xt- 

The above results have been obtained for mt = v^m^, the value employed in Ref. P]. 
In Table I, we compare the calculation of Ap using the formulation of Ref. with our 
own resonance and G.F. evaluations, over the range 120 GeV < m* < 220 GeV. We 
have checked that the numbers given in the third and fourth columns of Table I do not 
change when we identify Vq in Eq. (|^) with the zero of the integrand, i.e., the point 
where the resummation of {'KCpag/v)'^ terms matches the perturbative expression. It 
is apparent from Table I that the general features described above hold over the large 
range 120 GeV < mt < 220 GeV. In fact, the three calculations are similar in magnitude. 
Moreover, the formulation of Ref. gives results somewhat larger but quite close to our 
resonance calculation, the agreement being particularly good at low mt values, where the 
resonance picture is expected to work best. 

In order to illustrate the stability of the results with respect to a change of the scale, 
Pperh employed in the perturbative 0{as) calculations, in Table II we show the values of 
5{Ap)a^ (cf. Eq. (|ll])), 5{Ap)thr, and their sum for m^ = y/Sm^z and ppert = mt/^, m,t, 2mt. 
Here 6{Ap)thr is calculated on the basis of Eqs. (^ and (|), with as(mj) replaced by 
asifJ^pert) in the last term of Eq. (|^), and Vq chosen as the zero of the integrand (the factor 



(1 — t'^/3)/(l — v'^)'^ discussed before Eq. (||) is also appended). We see that there are 
variations in 5{/S.p)thr of —10% to 15%, which are not particularly large. Interestingly, 
they partly compensate the corresponding variations in (5(Ap)q^. Indeed, the overall QCD 
correction, ^(Ap)^^ + 6{Ap)thr, which is the physically relevant quantity, changes by only 
—3% to 4%, a remarkably small variation. 

3 Bohr-atom estimate and other observations 

It is instructive to make a simple estimate of the threshold effects on Ap by using an 
elementary Bohr-atom model of toponium |]l| . We have already employed this model to 
show that it leads to values of |-Rio(0)| within 20% of those obtained with the Richardson 



potential [|T6]. Now we want to apply it to illustrate the order of magnitude of S{Ap) 



thr- 



Then, instead of Eq. (H), we have 



5(ApW = -x,7rC|E^i%^, (12) 

n ^ 

where A;„ = Ciras(A;„)mi/(2n) is the momentum of the top quark in the nS orbital of 
the Bohr-atom model. We note that in this elementary estimate we have evaluated a^ 
at scale A;„. This is a simple generalization of Ref. |1[], where, for the ground state, a^ is 
evaluated at ki = Cpas{ki)'mt/2. In particular, for rrit = y/Smz and as{mz) = 0.118, 
we have ki = 16.7 GeV and as{ki) = 0.159. In the rough estimate of Eq. (p!2D, we have 
also disregarded the continuum enhancement above threshold. Because fc„ scales as l/n, 
as{kn) increases with n. We have iteratively evaluated kn and as(^n) up to ra = 50 and 
found Y^^=i(y^{kn)/n^ = 5.56 x 10~^, which is 1.38 times the IS" contribution. The sum 
converges rapidly; the first 12 terms already yield a factor of 1.36. Taking this to be an 
approximate estimate of the enhancement factor due to the resonances with n > 2, and 
normalizing Eq. (|l^) relative to Eq. (0), we have 

5{Ap)thr _ 2i3 Q^«(^i) (13) 

S{Ap)a, ' as{mt)' 

For rrit = y/Smz and as{mz) = 0.118, we have as{mt) = 0.109, ^^(/ci) = 0.159, and 
Eq. ([T3|) gives 42%, which is rather close to the results obtained from the detailed reso- 
nance approaches, i.e., Eqs. (|) ,(PDD divided by Eq. ([TT]). Although Eq. (^) is a rough 
estimate, it allows us to understand why the threshold effects on the p parameter, al- 
though nominally of 0{a^), can be as large as ~ 40% of the 0{as) contribution. Two 
factors are apparent: one is a large numerical coefficient, ~ 11, and the other is that 
the natural scale in the threshold contribution is ki <^ rrit, so that as{ki) is considerably 
larger than as{rnt). 

A relevant question is how these effects compare with unknown Ola"^) contributions. 
The leading 0{as) corrections to the p parameter are ~ 10% (cf. Eq. (pT])). If the same 
ratio holds between 0(a^) and 0{as), the threshold corrections we have discussed would 
be roughly 3 to 4 times larger for rrit ~ 160 GeV. It is known that the use of the run- 
ning top-quark mass absorbs most of the 0{as) corrections proportional to m^ in the p 



parameter and Z ^ bb amplitudes. If this was a general feature of the perturbative ex- 
pansion, one could control the bulk of such contributions and the threshold effects would 
neatly stand out. However, it is impossible to ascertain these features without detailed 
0{a1) calculations, which are not available at present. The effect of these threshold cor- 
rections on the electroweak parameters are not particularly large for rrit < 220 GeV. For 
example, for mn = 250 GeV and rrit = (130, 160, 200) GeV we found ||16[ shifts Aniw = 
-(42, 55, 77) MeV from the 0{as) contributions, Amw = -(14, 19, 27) MeV from thresh- 
old effects evaluated with the resonance method, and Amw = —(10,16,25) MeV from 
threshold effects evaluated in the G.F. approach. The ratio of threshold to 0{as) shifts 
in mw are similar but not identical to those we encountered in Ap. The reason is that 
the range of rrit values of interest is not really in the asymptotic regime, and the sub- 
leading 0{as) and threshold contributions to Ar |T^ (the relevant correction to calculate 
mw) have somewhat different rrit dependences. For example, for rrit = 160 GeV and 
rrifj = 250 GeV, the ratios for mw shifts are 0.35 in the resonance approach and 0.29 in 



the G.F. method |jT6[. For rrit = f^z (220 GeV), the corresponding ratios amount to 0.25 
(0.34) in the resonance approach and to 0.19 (0.34) in the G.F. method, with the latter 
giving somewhat lower values below rrit = 220 GeV. 

4 Conclusions and discussion 

1. We have pointed out that the conclusions of Ref. concerning the magnitude of 
the threshold effects discussed in Ref. |^ are without foundation. They are based on 
the consideration of a vector amplitude n(s) — 11(0), which, although very important 
in other applications, has a very small effect in the analysis of Ref. ||^. This should 
be obvious because heavy particles of mass m^ 3> s decouple in this amplitude. 

2. We have applied the analysis of '^mllthr{s') given in Ref. @] to study the threshold 
corrections to the p parameter for rrit = v^fnz and, after re-evaluating the quantities 
involved, found a result that is similar in magnitude to those we have reported 
previously. Contrary to the conclusions of Ref. [0, it is somewhat larger than 
our resonance calculation and significantly larger than our Green- function (G.F.) 
evaluation. In particular, when hard-gluon contributions are included, the result 
derived in the approach of Ref. is quite close to our own resonance calculation, 
well within the theoretical errors quoted previously P|. As illustrated in Table I, 
similar conclusions apply over the large range 120 GeV <mt< 220 GeV. We have 
also pointed out that the formulation of Ref. 0] leads, for rrit = v^^^z? to values 
of 5{Ap)a^ + 5{Ap)thr which are remarkably stable with respect to changes in the 
scale employed in the perturbative calculation (see Table II). 

3. In our opinion, the fact that the formulation of Ref. |^, when applied to Ap, gives 
results similar to ours (and, in fact, quite close to our own resonance calculations), 
supports the notion that these threshold effects can be reasonably estimated. By the 
same token, it does not support recent claims of large ambiguities in the threshold 



calculations 18 . 



4. It is worthwhile to point out that, when the corrected values of Spoie and 5thr obtained 
in the present paper are applied, the threshold effects in n(s) — 11(0) amount to 
~ 24% of the perturbative 0{as) corrections evaluated at scale rrit. Here rrit = 
\/3mz and as{rnz) = 0.115 have been assumed and, for simplicity, the hard-gluon 
correction is not included. (The numerical evaluation reported in Ref. ^ gives, 
instead, 15%). These numbers are somewhat smaller than the effects we encountered 
in the Ap case: 34% in the G.F. approach and 42% in the resonance framework (cf. 
Eqs. (^), (0) divided by Eq. (pJ])). However, we find it neither extraordinary nor 
unusual that radiative-correction effects may vary by factors of 1.42, 1.75 or, for 
that matter, 2.8, when applied to very different amplitudes. In particular, the logic 
behind the conclusion of Ref. seems rather strange to us. It is apparently based 
on the curious argument that a 15% correction in H(s) — H(0) is considered to be 
very small and that, as a consequence, a 34% or 42% effect in Ap is intolerably 
large. 

5. It is well known that, for large rrit, the widths of the individual top quarks become 
larger than the 15-25' mass difference, so that the bound-state resonances lose 
their separate identities and smear into a broad threshold enhancement 0. For 
this reason, we have expressed a preference to use the resonance formulation for 
rrit < 130 GeV and the G.F. approach for rrit > 130 GeV 0. However, we do 
not regard either method as "rigorous" and, in fact, in Ref. |^ we have assigned 
an estimated 30% uncertainty to their evaluation. In the analysis of electroweak 
parameters, we have found that the resonance approach of Ref. Q and the G.F. 
formulation of Ref. lead to similar results over a wide rut range, with the latter 
giving somewhat smaller values for rrit ^ 220 GeV. 

6. In order to make the relative size of the threshold effects more readily understand- 
able, we have presented a simple estimate based on an elementary Bohr-atom model 
of toponium (cf. Eq. (|13|)), and briefly commented on the possible magnitude of 
0{a1) corrections. 

7. We have not discussed here the non-perturbative contributions connected with the 
existence of a gluon condensate, since they are known to be exceedingly small in 



the case of the tt threshold [IS]; this has also been noticed in Ref. 
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